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Abstract — We consider a multi-hop wireless sensor network 
that measures sparse events and propose a simple forwarding 
protocol based on Compressed Sensing (CS) which does not 
need any sophisticated Media Access Control (MAC) scheduling, 
neither a routing protocol, thereby making significant overhead 
and energy savings. By means of flooding, multiple packets with 
different superimposed measurements are received simultane- 
ously at any node. Thanks to our protocol, each node is able 
to recover each measurement and forward it while avoiding 
cycles. Numerical results show that our protocol achieves close 
to zero reconstruction errors at the sink, while greatly reducing 
overhead. This initial research reveals a new and promising 
approach to protocol design through CS for wireless mesh and 
sensor networks. 

I. Introduction 

Wireless sensor networks are used in large sensing areas 
where multi-hop routing is needed to forward the information 
from the sensors to the sinks. When a sensor wakes up 
according to an event and then measures a physical value, 
it first needs to schedule its transmission in an interference- 
free slot using a Media Access Control (MAC) layer scheme 
and then to use a specific routing protocol in order to fix the 
next hop before reaching the final sink. Those steps include 
heavy operations for small sensors and consume a big quantity 
of energy. However, they are mandatory in order to avoid 
collisions and to determine the best routes to reach sinks. 

Recently, the groundbreaking theory of Compressed Sensing 
(CS) was developed, which stipulates that a vector with 
correlated entries, i.e., that can be transformed into a sparse 
vector through a transformation basis, can be recovered with 
high accuracy from a few random projections onto another, 
incoherent basis [1] [2] [3], CS, which has been widely used 
in the domains of digital signal and image processing, is 
envisioned as a highly promising tool for improving the per- 
formance of wireless communication networks. In particular, 
based on the decentralized compression techniques presented 
in [4], there has been a number of works exploiting CS in 
wireless sensor networks dealing with space-time correlated 
data as in [5] [6]. In these works, the underlying assumption is 



that all sensor nodes have a data to send which are gathered 
at the sink and reconstructed with few observations using 
CS techniques. The transmission of each data follows some 
predetermined routes, such that packets are received without 
errors. In [7], the CS technique is combined to the well-known 
Network Coding (NC), which provides an efficient method 
for data communication in wireless networks, as multiple 
packets aggregated in a single packet are decoded using prior 
information. Introducing CS enabled to mitigate the problem 
of exploding header size with NC where the header grew 
proportionally to the number of aggregated packets, since 
it became possible to decode the salient information from 
fewer number of packets. Note also that, as a consequence, 
CS enables to significantly decrease the delay required by 
NC for gathering packets, as at least n different combined 
packets were needed to decode n packets. However, in this 
work too, the transmissions follow predetermined routes in 
order to avoid any collision, i.e., packets are aggregated 
from node to node and not "over-the-air". By contrast, [8] 
considers a one-hop sensor network where multiple nodes may 
transmit their measurement simultaneously to the sinks. As the 
measured events are assumed to be sparse, i.e., the number 
of measurements simultaneously received at the sink is much 
smaller than the total number of sensors, the sink is able to 
recover each measurement from the superimposed signals from 
few observations, thanks to CS algorithms. Such an approach 
is also taken for developing new multiple access schemes for 
random data traffic in [9], for channel reservation in [10], and 
for downlink scheduling in [11]. 

In this work, we consider a multi-hop wireless sensor 
network and take the approach of [8] — [ 11] where the events 
to be reported occur sparsely. Each sensor with a value to 
report combines it with its signature sequence giving its 
Identity (ID) and broadcasts the resulting vector to its neighbor 
nodes. We develop a protocol based on CS and flooding that 
enables the sink to obtain and reconstruct the sparse mea- 
surement data with high accuracy without any heavy routing 
nor MAC protocol. In a one-hop network, the simultaneously 



received measurements may be resolved by the CS-based 
algorithms in [8]— [ 11], but they pose major problems in a 
multi-hop network where every node forwards all received 
packets, as the number of interfering measurement packets 
may be increased drastically, thereby causing the CS-based 
algorithms to perform poorly due to loss of sparsity. Moreover, 
if different measurements of a same source but generated 
at different times are contained in the same packet, they 
would be hardly resolvable. To alleviate the aforementioned 
issues, in our protocol, the length of the signature sequences is 
further increased so that each sequence can identify the source 
node and time stamp, which allows to avoid cycles in the 
flooding process as follows. When neighboring nodes receive 
the superimposed packets, they first recover each measurement 
data from an under-determined system of equations, using a 
CS algorithm based on £1-^2 optimization, commonly used 
in the case of noisy measurements [12]. After decoding each 
data, they check in their local tables if it was already sent, 
using the recovered origin node ID and time stamp from 
its sequence; if it is the case, they discard the data. If not, 
they again superimpose each data into a unique packet which 
is broadcasted to their neighbors. This process is repeated 
until reaching the sink. In addition to avoiding cycles, this 
process allows to maintain the "sparsity" of the superimposed 
data by discarding superfluous packets, thereby guaranteeing 
a good performance of the CS algorithms. Still, the length 
of these sequences is far below the total number of possible 
combinations of IDs and time stamps, which would be required 
in a CDMA-based protocol. 

Many advantages are offered by our new protocol: 

• Sensors do not need any MAC scheduling that consumes 
resources, energy and increase the delivery delays. 

• Sensors do not need any routing protocols that need 
many control messages that waste the bandwidth and also 
increase the energy consumption. 

• Decoding in the sensors needs processing power and 
consumes energy. However, those operations represent a 
negligible consumption compared to the radio reception, 
listening, transmission used by the MAC scheduling. 

• The time stamp represents an additional overhead but 
is bounded by, e.g., the maximum number of hops in 
a network. We show that it represents a negligible part 
compared to the overhead induced by a routing protocol 
and even compared to a flooding-based protocol with 
CDMA. 

• While discarding packets that have already been for- 
warded for keeping sparsity of the superimposed data 
and avoiding cycles, our protocol still offers diversity for 
each measurement data as a packet may reach the sink 
by many routes. Even if a packet is lost, the sink or the 
intermediate nodes may receive it by other routes. 

The paper is organized as follows: after giving some basic 
results in Compressed Sensing and presenting related works, 
we explain the system model and our proposed protocol. The 
simulation results show its effectiveness, while reducing the 



amount of overhead compared to conventional protocols. 

II. Results from Compressed Sensing 

Compressed Sensing is a newly developed mathematical 
theory that enables to solve underdetermined systems of 
equations under the sparsity prior of the solution, a problem 
formulated as 

(Po) : min||x||o, subject to y = Ax, (1) 

X 

where x e K w , y e R M , A £ R MxN with M < N, and 
||x|| = card{i,Xi ^ 0}, where sparsity of x is defined by 
||x|| = K « N. However, this is a non-convex optimization 
problem that requires combinatorial search over all possible 
sparse combinations of columns of A, which rapidly becomes 
intractable. Therefore, convex relaxation of this problem has 
been considered by the following £1 -minimization, 

(Pi): min||x||i, subject to y = Ax, (2) 

X 

where ||x||i = J27=i \ x i\> f° r which various efficient algo- 
rithms have been developed, among which Pursuit algorithms 
and Iterative Shrinkage Algorithms [12]. 

Under the assumption of noisy observations, the optimiza- 
tion problem may be reformulated as 

(Pi): min||x||i, subject to ||y — Ax|| 2 < e, (3) 

X 

which can be equivalently expressed as an ii —£ 2 minimization 
problem, 

(P*_ 2 ): minA||x|| 1 + i|| y -Ax||i, (4) 

x I 

where the parameter A can be viewed as a Lagrange Mul- 
tiplier making a trade-off between representation error and 
sparsity of the solution. It was shown that random matrices 
A whose entries are drawn from Gaussian or Sub-Gaussian 
(e.g., Bernouilli) distributions, guarantee stable recovery in the 
noisy case, provided that the number of measurements obey 

M >cKlog(N/K), (5) 

where c is a constant (this condition guarantees that A 
will guarantee the Restricted Isometry Property (RIP) with 
high probability, enabling stable recovery by ii -minimization, 
see [12] and references therein). In our proposed protocol, we 
will employ the Iterative Shrinkage-Thresholding Algorithm 
(ISTA) [13] as it efficiently solves the li — £2 minimization 
(P 1 A _ 2 ) with very low complexity. 

III. Related work 

In the multiple access schemes of [10] [9] [11] for a single- 
hop system, several Mobile Stations (MS) access the Base 
Station (BS) at the same time, on the same frequency. Assum- 
ing J MSs in total, each MS j is assigned a pseudo-random 
signature sequence aj, a vector of size M, which entries are 
generated by Bernoulli random variables of probability 1/2, 
with values ±1. Thus, in the absence of noise, the received 
signal y at the BS may be expressed as 

y = [ai ... aj ... a, 7 ] v = Av, 



where v is the vector of size J with component {v}j = 1 
if MS j transmitted, and {v}j = if not. As the number of 
transmitting MSs is assumed to be much smaller than the total 
number of MSs J, v is a sparse vector. Hence, v is recovered 
at the BS using CS algorithms solving problems (0, or in 
the noisy case. Note that, the usage of these pseudo-random 
signature sequences to identify each MS may be regarded 
as similar to a Code Division Multiple Access (CDMA) 
system. However, CDMA systems typically require orthogonal 
codes to differentiate each MS, which would require each 
sequence &j to be at least of size J, whereas these multi- 
access systems with CS deal with non-orthogonal sequences 
since M « J, which can be viewed as an overloaded CDMA 
system, in which case traditional CDMA receivers such as 
the linear minimum mean squared error filter perform very 
badly, as shown in [11]. Nevertheless, we will also compare 
our proposed protocol to a CDMA-based reference protocol 
that does not require any routing. 

IV. System Model 

We consider a multi-hop wireless sensor network with N 
nodes and one sink, forming a lattice as shown in Fig. Q] Each 
sensor node S n can be either in transmit or receive mode 
at any one time. If a sensor detects an event, the network 
forwards the corresponding measurement to the sink based on 
the proposed protocol. As in [8] which considers a one-hop 
sensor network (but multiple sinks), we assume that there may 
be up to K events occurring simultaneously within the whole 
network, but these events are considered to be sparse compared 
to the number of nodes, namely K « N . 

0^0-0-0 
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Fig. 1. Multi-Hop Wireless Sensor Network Model, Step 1 

We assume digital transmissions as in [8], i.e., measure- 
ments x E { — 1,1}, instead of directly transmitting scalar 
values by amplitude modulation as in [4], since interfering 
analog signals pose difficult problems in a multi-hop network. 
For example, a node could simultaneously receive two packets 
from different paths, but containing the same measurement x, 
in which case it would erroneously decode 2x. As the goal 
here is to present the benefits of our protocol, we consider the 
digital case for simplicity, but the protocol will be developed 
to accommodate analog transmissions in the next phase. 



Given the lattice structure of the network, it is assumed 
that each node only communicates with its closest neighbors 
at distance d, so there may be 2, 3 or 4 neighbors depending 
on the node's location. For example, a packet sent by node Si 
in Fig. Q]will be received by nodes S% and S§ only, while a 
packet sent by node SV will be received by nodes S2, Sq, 
Ss and S12. For simplicity, all channels are considered to 
be Additive White Gaussian Noise (AWGN), and all links 
between neighbors at distance d have the same gain. Channel 
fading and path loss effects will be considered in the future 
worlo Note that all nodes are synchronized as in [4] [8], but 
there are interesting design issues in the asynchronous case, 
to which the protocol may be extended. 

P2 

0-0^0 
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0^0^0 
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Fig. 2. Multi-Hop Wireless Sensor Network Model. Step 2 



V. Proposed Protocol 

The proposed routing protocol consists in a simple flooding 
whereby upon packet reception, each node broadcasts this 
packet locally and then all its neighbors forward this informa- 
tion to their neighbors and so on, until it reaches the sink. In 
this setting, each packet will be forwarded through different 
paths, creating multiple copies that provide diversity at the 
sink. As in the works presented in Section [TTU we use in 
our system the Bernoulli sequences for identification, forming 
the independent projection vectors required for recovering the 
compressed signals. However, if each source node is only 
assigned a unique signature sequence as in the previous works, 
it will not be possible to distinguish whether the multiple 
received packets are due to self-interference, i.e., they are 
copies of a same original packet, or whether they correspond 
to measurements from the same source but taken at different 
times. 

To alleviate this problem, we associate each measurement 
with a unique sequence a n ,i, vector of size M, that identifies 
the pair (origin node ID S n , number of hops I), where S n 
is the sensor that originally detected this measurement, and I 
is its time stamp parameter, defined as the number of hops 
elapsed by forwarding the packet bearing this measurement, 

'Note that, although fading and path loss effects are considered in [8], 
the channel gains between each sensor and sink are assumed to be known 
perfectly as they form the columns of sensing matrix A. 



counted from its origin node up to its current location in the 
network. We have a n j G { — 1,+1} M , given by a pseudo- 
random noise signature sequence generated by Bernoulli ran- 
dom variables of probability 1/2, with ±1 values. Thus, fixing 
a maximum delay for a packet or time to live L in number 
of hops, which depends on the network size, there are NL 
sequences a n j of size M, where we set M < NL for 
reducing overhead consumption. These sequences are assumed 
to be known at each node. They are gathered into matrix 
A of size M x NL, whose column q = (n — 1)L + I 
contains a„ /. For notational convenience, let x denote the 
iVX-lengfh sparse vector where {x}( n ~i)L+i = x k,i which 
groups all the measurements with possible hop counts, i.e., 
x = [xi.o---Xi,L---Xn,i---XNft---%N.L\ T where the L + 1 first 
terms store the measurements of source node S\ for I = 0..L 
and similarly, the L + 1 last terms for source node Sn for 
l = 0..L. 

Then, the proposed algorithm works as follows: 

1. If source node Sj in transmit mode detects a new event 
with measurement value Xj,o 6 { — 1, +1} initialized by / = 0, 
he transmits the packet of M bits 

Pi = x jfl*j,o, (6) 

which is received by all its neighbor nodes k € N(Sj). For 
example, assuming Si and S3 to be source nodes in Fig. Q] 
packets pi = £1,081,0 and p% = £3,083,0 are received by their 
respective neighbors, i.e. S2,Sq for Si and S2, S4, Sg for S3. 

2. At a time t, if sensor node Sk is in receive mode, he receives 
signal yk containing all simultaneous transmissions from his 
neighbor nodes n G N(Sk)- Thus, the received signal at a node 
Sk is written 

Yfc = ^2 Pn=^2 x J,l a jd+ z k, (7) 

ri£N(S fc ) j,l 

where Xj,i, the measurement with source node Sj is included 
in p n if it was actually forwarded by node S n . Vector Zk 
of size M, denotes the AWGN. Using the matrix A defined 
above, this is equivalently reformulated as 

y fe =Ax fc +z fc , (8) 

where X& is the vector of size NL defined as 

, , _ J {x}, ; if {x}i = xj.i is received by Sk 
l x fcj 4 - I o otherwise 

In the example of Fig. Q] note that Sj = S n , i.e., the source 
nodes and transmit nodes coincide since we describe Step 1, 
so the received signal at receiver node 5 2 is written y 2 = 
Ax 2 + z 2 , where x 2 = [xi t o 0...0 2:3,0 0...0] T . 

3. Receiver node Sk decodes each superimposed measure- 
ments contained in the received signal y^, using CS princi- 
ples. In particular, this problem can be formulated by I1-I2 
optimization as in 

x fe =axgmin-||y fe -Ax fe ||| + A||x fc [|i, (9) 

which we solve by ISTA with very good reconstruction 
abilities with low complexity. 



4. The reconstructed vector Xfe is renormalized so that each 
component belongs to { — 1, +1}. 

5. After decoding the received measurements at node Sk, 
their respective (source ID, time stamp)-sequence &jj could be 
identified. Node Sk then compares the sequence aj.i of each 
newly decoded measurement to the sequences of previously 
received ones that are stocked in its local table Q(Sk) = 
{a qt i> received by S k }- 

• if there exists a sequence a g .c in Q(Sk) with q = j and 
/' < I: the same data was already received before with 
high probability, since no fading is assumed, such that 
any new measurement should have a hop count / < 
unless the packet that followed the shortest path was lost. 
However, this will rarely occur in this setting. Thus, Sk 
discards the decoded measurement Xj^ajj, from the set 
of data to be forwarded next. 

« if there exists a sequence a q j/ in Q(Sk) with q = j and 
I' > I: although both share the same origin node Sj, this 
is a new data as it has a smaller hop count. Therefore, 
Sk will forward the decoded data Xjj&jj. 

• in all other cases, Sk will forward Xj.iSjj. 

6. For all data Xj^Sj^ to be forwarded, the hop count is 
incremented to I + 1. All decoded measurements as well as 
sequences in Q(Sk) for which I > L are discarded. Then, Sk 
superimposes all the remaining data to be forwarded into a 
packet pk composed of M bits, 

p k = Ax fc = y^x 3 ya.,y, (10) 

3,1 

and forwards it. Thus, in Fig. [2] node 5 2 sends p 2 = ii,iai,i + 
£34834. 

7. The sink node runs the reconstruction algorithm for all in- 
coming packets during the session of duration T out . Whenever 
multiple versions of the same data are received, more diversity 
gain may be achieved to improve the decision accuracy. 

VI. Numerical Results 

We first evaluate our proposed protocol in a network of 25 
nodes. In a session, K sources generate measurements that are 
forwarded to the sink as explained in Section [V] These sparse 
generation events occur at random times during each session, 
e.g., they may occur simultaneously. Each session ends after 
time out T out in hop counts, at which point the measurements 
collected by the sink are evaluated. The Signal-to-Noise Ratio 
(SNR) between neighboring nodes is fixed to 10 dB. Each 
simulated point is averaged over 300 sessions. We evaluate 
the reconstruction error as in [5] [6], averaged over all sessions, 
defined for one session as 



where Xo is the vector of size containing the original mea- 
surements for all N sensors, i.e., its corresponding components 
are in { — 1 , + 1 } if sensor Sj is an origin node, and zero 
otherwise, while x denotes the reshaped vector of collected 
measurements at the sink at the end of a session. 
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Fig. 5. Proposed Protocol: Reconstruction Error performance for different 
values of A and s = 6, 8, 10; T out = 30, N = 49, M = 40 and L = 7. 
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Fig. 4. Proposed Protocol: Reconstruction Error performance for different 
values of M, A and Tout = 30, N = 49, L = 7. 



Fig.[3]shows the average reconstruction error of our protocol 
for K — 2, 10, for sequence length M — 30, L — 5 and 
Tout = 20, 30. We fix A = 0.5 in the decoding algorithm by 
ISTA. We observe that in both cases, excellent reconstruction 
performance is achieved for all numbers of event sources K, 
even for larger values of K which imply lower sparsity in 
the number of superimposed measurement data. The error 
performance is improved by increasing T out , as more diversity 
gain is achieved at the expense of delay. 

Next, we consider a network of N = 49 nodes for T ou t = 
30, L — 7 and M = 40,50. For fixed A = 0.5, the 
reconstruction error increases as K increases, due to lower 
sparsity: as the number of interfered packets increases, it 
becomes more difficult to correctly receive and recover all 
measurements at each node. However, by increasing M from 
40 to 50, the error performance decreases significantly for 
K = 8, 10, since it provides a larger number of independent 
observations that allow to correctly decode a higher number 
of superimposed data by improving the estimation accuracy 
of the CS algorithm. Note that, even for small values of K, 
the number of interfering packets to be decoded at a certain 
node may be much more, since multiple copies of a given 



packet may be received through different paths, reducing the 
sparsity of the received signal y>. However, increasing M 
creates a larger overhead, even though it is far smaller than 
what is required by conventional protocols, as shown later. 
Instead of increasing M, Fig. H shows that the reconstruction 
error can be made close to zero by adapting the parameter A to 
A op i, given the number of event sources, namely the sparsity 
of the system. Thus, Fig. [5] evaluates the reconstruction error 
performance for the number of event sources s = 6, 8, 10, 
and varying values of A. Note that curves for the "sparser" 
cases s = 2,4 are not represented, as their reconstruction 
error stays close to zero for A £ [0.4, 0.9]. We observe 
that, as the number of sources increases, and hence as the 
sparsity diminishes, the reconstruction performance by ISTA 
becomes more sensitive to an accurate optimization of A. 
In this case, the optimized values of A could be chosen as 
A op t = 0.6, 0.7, 0.8 for s = 6, 8, 10, respectively. The tendency 
for X op t to increase with s can be interpreted by the fact that 
more weight should be put on the l\ term in (|4]i for finding a 
sparse solution, as more and more "denser" candidate solutions 
appear. Of course, the excellent reconstruction performance 
achieved by X opt comes at the price of increased computational 
complexity at each node, even though it should remain at a 
reasonable level given the very low computational complexity 
of ISTA. On the other hand, the alternative to optimizing 
A is to sufficiently increase M, which also provides very 
good performance as already shown in Fig. [4] with M = 50, 
A = 0.5. The trade-off is thus between algorithm complexity 
or increased overhead, which is evaluated next. 

We compare the proposed protocol with two reference 
algorithms: a Conventional (Conv.) algorithm based on AODV 
routing [14] and MAC scheduling, and a flooding-based pro- 
tocol similar to our proposed algorithm, but using CDMA 
as explained in Section [Til] The comparison is only made 
in terms of overhead, since both reference algorithms should 
ideally provide excellent reconstruction performances under 
this system model, given the absence of interference for the 
first one, and due to the orthogonal codes for the second. 
As overhead and data (origin node ID and measurement) are 



not separable in the proposed protocol, we evaluate the total 
amount of packets generated by each method. For the Conv. 
algorithm, the bounds are roughly given by 



D 



AODV — B RRE q(N — 1) + B RREP h, 



(12) 



where B RR eq = 32 bytes for Route Request, B RR ep = 28 
bytes for Route Reply, and 



B 



MAC 



'ID+DATA 



?CRC) 



2.5/1, 



(13) 



where bio+DATA = 11 bits and bene — 8 bits for Cyclic 
Redundancy Check. The lower bound B™™ v is given by h = 1 
with the origin node at 1-hop distance from the sink, and the 

upper bound B™^, b Y h = 4 for N = 25 ( h = 6 for N = 49 X 
with the origin node at a corner of the lattice network, and 



Tjmin r>min , r>min 

n Conv — D AODV "+" n MAC> 

r>max rymax , r>max 

n Conv — & AODV "+" n MACi 



(14) 
(15) 



for each source measurement. For the proposed protocol, the 
number of bits obtained by simulations, may be expressed by 
Bprop = MPtotatProp, where Ptotai,Pro P is the total number 
of packets occuring in a session, including all the forwarded 
ones due to flooding. Similarly, for the reference scheme with 
CDMA, we have B C dma = N LP to tai,CDMA, as each se- 
quence should have a minimal length of NL to guarantee near- 
zero reconstruction error as explained in Section [HI] Fig. [6] 
shows that the proposed protocol largely decreases the amount 
of packets per session in [kbytes] for networks with 25 and 
49 nodes, corresponding to 80 to 85% reduction compared to 
the conventional algorithm with routing and MAC. Compared 
to the reference scheme with CDMA, although the gain of our 
protocol is around 75% for N = 25, it becomes tremendously 
high for N = 49, going up to 90%. Due to the large amount of 
bits required for having orthogonal sequences in the CDMA 
system, the total amount of packets drastically increases and 
even surpasses the conventional routing and MAC based 
protocol as the network enlarges, showing evident scalability 
issues that are alleviated by our proposed protocol. Thus, the 
proposed protocol provides very low reconstruction errors, 
while achieving tremendous overhead savings that become 
even higher as the network enlarges. 

VII. Conclusion 

We have proposed a novel protocol enabling non-scheduled, 
simultaneous transmissions for delivering measurement data 
to the sink in multi-hop wireless sensor networks. Due to 
the sparse nature of the sensed events, the superimposed 
measurements resulting from interfered packet reception at 
each node can be recovered by CS algorithms. By making use 
of time stamps and signature sequences for each measurement, 
our protocol achieves low reconstruction errors, bringing sig- 
nificant savings in terms of overhead and energy, as compared 
to conventional routing and MAC scheduling strategies, as well 
as CDMA and flooding based protocol. 

This fundamental work opens new vistas for further re- 
search. The protocol will be extended to more general sys- 
tems, including path loss and channel fading, various network 
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Fig. 6. Average amount of transmitted packets per session in [kbytes] for 
two cases: N = 25, M = 30, L = 5, T out = 30 and N = 49, M = 50, 
L = 7, T out = 30. 

topologies, asynchronous nodes and, the more challenging 
case of analog transmissions. Another promising direction is 
the design of new algorithms based on the proposed ideas 
for general multi-hop mesh networks in the context of uni- 
cast/multicast transmissions. 
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